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Abstract 



Bose-Einstein Correlations in one and two dimensions have been studied, with 
high statistics, in charged current muon-neutrino interaction events collected with 
the NOMAD detector at CERN. In one dimension the Bose-Einstein effect has 
been analyzed with the Goldhaber and the Kopylov-Podgoretskii phenomenologi- 
cal parametrizations. The Goldhaber parametrization gives the radius of the pion 
emission region Rq = 1.01 ± 0.05(stat)loQg(sys) fm and for the chaoticity pa- 
rameter the value A = 0.40 ± 0.03(stat)^QQg(sys). Using the Kopylov-Podgoretskii 
parametrization yields Rrp = 2.07 ± 0.04(stat)+Q°4(sys) fm and Xj^p = 0.29 ± 
0.06(stat)^QQ4(sys). Different parametrizations of the long-range correlations have 
been also studied. The two-dimensional shape of the source has been investigated 
in the longitudinal co-moving frame. A significant difference between the transverse 
and the longitudinal dimensions is observed. The high statistics of the collected 
sample allowed the study of the Bose-Einstein correlations as a function of rapid- 
ity, charged particle multiplicity and hadronic energy. A weak dependence of both 
radius and chaoticity on multiplicity and hadronic energy is found. 



Deceased 
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1 Introduction 



The quantum mechanical wave function of two identical bosons has to be symmetric under 
particle exchange. The symmetrization gives rise to an observable interference pattern 
which enhances the number of identical bosons emitted close to one another in phase space. 
Such Bose-Einstein Correlations (BEC) were observed for the first time in astronomical 
measurements of photon pairs emitted by stars [1] and soon after for like-sign hadrons 
produced in pp annihilations [2]. Since then BEC, were also measured in several other 
types of particle interactions (for a review see [3]). The shape of the BEC depends on the 
spatial and temporal distributions of the boson source and on its degree of coherence. The 
theoretical aspects of the BEC were developed in the papers of Kopylov and Podgoretskii 
[4] and Cocconi [5]. From these studies it appears that the measurements of BEC may be 
important to gain an understanding on the dynamics of the particle interactions yielding 
like-sign bosons in the final state. 

Previous measurements of the BEC effects in neutrino interactions have been performed 
by the Big European Bubble Chamber Collaboration (BEBC) [6] including data collected 
on a variety of targets by both BEBC at CERN and the 15-foot Bubble Chamber at 
Fermilab. Nevertheless, the number of events globally collected by these experiments is 
still about one order of magnitude smaller than the data set collected by NOMAD and 
used in this paper. 



2 The phenomenology of BEC 

The BEC effect can be parametrized in terms of the two particle correlation function R 
defined as: 



where pi^ are the particle four- momenta, D(jpi,p2) is the measured two-particle density 
and Dq(pi,P2), the particle density in the absence of BEC. Do(pi,p2) should include 
any other two-particle correlations such as those coming from phase space, long-range 
correlations, charge effects, etc. which in the ratio should be divided out leaving only 
the BEC effects. According to the Goldhaber parametrization [2], which assumes that 
the emitting sources of identical bosons are described by a spherical Gaussian density 
function, BEC are usually parametrized as: 



where Q 2 = —(P1—P2) 2 = M^—Am 2 , with M V7T the invariant mass of the pion pair, Rq the 
width of the Gaussian distributed emitting source and the pion mass. The chaoticity 
(or incoherence) parameter A measures the degree of coherence in pion production, i.e. 
the fraction of pairs of identical particles that undergo interference (0 ^ A ^ 1). 



R(puP2) = D(p 1 ,p 2 )/D (pi,p 2 ) 



(1) 




(2) 
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The Kopylov-Podgoretskii (KP) parametrization [4] corresponds to a radiating spherical 
surface of radius Rrp with pointlike oscillators of lifetime r: 



R(Qt,Qo) = 1 + A AJl{Q t R KP )/{Q t R KP ) 2 I 1 + (Q r) 



(3) 



where J\ is the first-order Bessel function, p = pi + p 2 , Q = p\ — p 2 , Qo = \Ei — E 2 \, 
Qt = Q x P I \v\ ■ This parametrization is not Lorentz invariant and the variables are 
calculated in the centre of mass of the hadronic final state. It can be shown that at small 
values of Q t and Q the relation R KP xs 2R G is expected [6]. 



The shape of the hadronic source can be measured by studying BEC as a function of 
the components of the vector Q. It is convenient to perform this study in the so-called 
longitudinal centre of mass system (LCMS). This reference system is defined for every 
particle pair as that where p = p\ +p 2 is perpendicular to the axis defined by the hadronic 
jet direction (see fig. 1). With this choice, possible effects caused by the Lorentz boost are 

— * 

avoided. In the LCMS, Q is decomposed into the following components: Qi on g, parallel to 
the hadronic jet axis; Qt,out, collinear with p and the complementary Qt, S ide, perpendicular 
to both Qi on g and Qt,out- I n this analysis we use the longitudinal component Q\\ = Qi ong 
and the perpendicular component Q± = ^ Ql tOUt + Qt, S ide ( see fig-1)- The parametrization 
of the correlation is then performed separately for the longitudinal Q\\ and transverse Q± 
components as suggested in ref [2]: 

R(Q\\ ,Q±) = 1 + Aexp ( - Q\R\ - Q\R\) (4) 

the longitudinal and transverse dimensions of the hadron source being represented by Rn 
and R±, respectively. 



3 Experimental procedure 



3. 1 The NOMAD experiment 



The main goal of the NOMAD experiment [7] was the search for — > v T oscillations in 
a wide-band neutrino beam from the CERN SPS. The full data sample, corresponding 
to about 1.3 million charged- current (CC) interactions collected in four years of data 
taking (1995-1998) in the detector fiducial volume, is used in the present analysis. The 
data are compared to the results of a Monte Carlo simulation based on modified versions 
of the LEPTO 6.1 [8] and the JETSET 7.4 [9] generators for neutrino interactions and on 
a GEANT 3.21 [10] based program for the detector response. BEC effects are not included 
in the Monte Carlo. For the analysis reported below we have used a Monte Carlo sample 
of size comparable to the data. 
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Fig. 1. The definition of the LCMS system. 



3.2 The NOMAD detector 
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Fig. 2. NOMAD Apparatus. 



The tracking capabilities of the detector are essential for the study of BEC. The NOMAD 
detector shown in fig. 2 is especially well suited for this. It consists of an active target 
of 44 drift chambers, with a total fiducial mass of 2.7 tons located in a 0.4 T dipole 
magnetic field. The drift chambers (DC) [11], made of low Z material (mainly carbon), 
serve the dual role of a nearly isoscalar target for neutrino interactions and of a tracking 
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medium. These drift chambers provide an overall efficiency for charged particle recon- 
struction greater than 95% and a momentum resolution which can be parametrized as 
~$t = VT ® where the track length L is in metres and the track momentum \p\ in 

GeV/c. This amounts to a resolution ^ 3.5% for \p\ ^ 10 GeV/c. Reconstructed tracks are 
used to determine the event topology (the assignment of tracks to vertices) and to recon- 
struct the vertex position and the track parameters at each vertex (primary, secondary 
V0, etc. . . ). A transition radiation detector (TRD) [12] is used for electron identification. 
The pion rejection achieved for isolated tracks is 10 3 with a 90% electron identification 
efficiency. A lead-glass electromagnetic calorimeter (ECAL) [13] located downstream of 

the tracking region provides an energy resolution of 3.2%/y^fGeV] © 1% for electromag- 
netic showers and is essential to measure the total energy flow in neutrino interactions. 
In addition an iron absorber and a set of drift chambers located after the electromagnetic 
calorimeter are used for muon identification, providing a muon detection efficiency of 97% 
for momenta greater than 5 GeV/c. 

3.3 Event selection 

The identification of u^CC events requires the presence of a primary negative muon in the 
final state, i.e. a track segment in the muon detector matched to a track reconstructed 
in the drift chambers. The muon momentum and its transverse component (relative to 
the neutrino beam) are required to be greater than 5 GeV/c and 0.5 GeV/c, respectively. 
Preliminary cuts are applied to ensure good quality event reconstruction 2 : 

• Number of primary charged tracks (excluding the muon) N ch ^ 2; 

• Muon energy > 5 GeV and hadronic energy E hadrons > 5 GeV; 

• Hadronic invariant mass W ^ 2 GeV (to reject quasi-elastic events and baryon reso- 
nance production); 

• Event vertex within the fiducial region of the DC target. 

Tracks to be used for BEC are then selected using the following criteria: 

• Only primary tracks are selected, i.e. either belonging or pointing to the primary vertex; 
in the latter case, to avoid potential dangerous contamination from photon conversions, 
the track first hit must occur no further than 15 cm downstream of the vertex along 
the detector axis (the z axis); 

• A minimum momentum of the track is required: \ptrack\ ^ 100 MeV/c; 

• A minimum number of hits is used to build the track: Nh its ^ 12; 

• A good momentum resolution is required: ^ 6%, where A \p\ is the the uncertainty 
on the momentum of that track calculated by the reconstruction program; 

• The track should not be identified as an electron by the TRD and the ECAL; 

• The track should not be identified as a proton by the range- momentum correlation 
method (see ref. [15] for details). 

2 additional details on the analysis can be found in [14]. 
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The total number of events selected by these cuts is 398K. These events contain 544K 
(++), 143K ( ), and 852K (+-) pairs. 

After the event and track quality cuts, we have performed a preliminary analysis of the 
simulated events to assess the purity of the track identification. The tracks investigated are 
those obtained from the full Monte Carlo simulation. An appropriate algorithm allowed an 
association between the reconstructed tracks and those generated at the primary vertex. 
We have found that the positive and negative samples of particles used for BEC studies 
contain respectively ~ 61% of tt + and ~ 77% of tt~ . The relative contributions of various 
positive and negative particles in NOMAD entering the correlation plots are listed in 
table 1. The tracks labeled "not recognized" are those for which no association with a 



Particle 


Percentage 


7T+ 


61.2% 


V 


18.4% 


K+ 


5.6% 




0.1% 


e+ 


0.2% 


Not recognized 


14.5% 




77.2% 


V 


2.4% 


K~ 


7.8% 


e~ 


0.3% 


Not recognized 


12.3% 



Table 1 

Composition of the charged particle sample used in the analysis. 

generated primary track was found. These tracks are produced by secondary interactions 
or photon-conversions. As we can see from table 1 the negative tracks exhibit a better 
pion purity than the positive ones. Contaminations from electrons and positrons which are 
mainly present in the "not recognized" samples could be dangerous since these particles 
come from photon conversions and therefore they can populate the low Q region where 
BEC effects are expected. This problem, together with the effects caused by the proton 
and kaon contaminations, will be discussed in section 4. 

The Monte Carlo simulation was also used to verify the experimental resolution in the 
determination of the BEC parameters. BEC effects are expected to occur in any of the 
kinematical variables Q, Q t , Qo, Q\\, Q± at small values of these parameters 0.2 GeV) 
and we have verified that the resolution in any of these variables is ^ 0.02 GeV. The BEC 
parameters are obtained by fits to the experimental distributions in the interval 0.0-1.5 
GeV. This interval is large enough to study possible long-range correlations as well. 
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4 The reference samples 



We have studied several alternatives for the choice of the reference sample D (pi,p2) 
used in eq. 1. In principle the Monte Carlo events, which do not contain BEC, would be 
good candidates. However the capability of the Monte Carlo to accurately reproduce the 
data (except for BEC) especially in the tiny phase-space region where BEC are present is 
limited and other methods based on the data themselves must be found. Several methods 
have been used in previous experiments to build the reference sample from the data (see 
for example [6]). They are: 

• The reference sample is formed of all unlike-sign pairs; 

• The reference sample is formed by building a so called "mixed event": a hadron from 
one event is combined with a hadron of the same charge, chosen at random from an 
another event that has approximately the same kinematical characteristics: total hadron 
momentum, hadron energy, charged multiplicity... 

• The reference sample is formed by pairing unlike-charge hadrons from the same event 
after the transverse momenta pt (with respect to the current direction) have been in- 
terchanged at random in the hadronic centre-of-mass system (c.m.s). 

We have carefully tested the three methods with a full Monte Carlo simulation (discussed 
in more detail later in this section) which includes also the response of the detector. The 
Monte Carlo results reproduce correctly the inclusive particle distribution in neutrino 
interactions, but not the correlations among particles. This Monte Carlo is therefore 
adequate to study the bias introduced by the reference sample. In fact, in the absence 
of BEC, the distributions in R should be flat or, in any case, have no structure at small 
Q (^ 0.2 GeV) which would distort the study of R. We have found that none of the 
three reference samples completely fulfills this requirement, however the unlike-sign one 
was eventually found to be the most adequate in the BEC region. Moreover this reference 
sample has been used by the great majority of previous BEC studies. 

BEC effects are then investigated by looking in the data at the following ratio: 

= "like-sign" pion-pairs = N ++ (Q) + iV--(Q) 
"unlike-sign" pion-pairs N+_(Q) 

The Monte Carlo samples are used to estimate possible spurious BEC effects from non- 
pion contaminations present in the sample as discussed in the following section. 



4-1 Systematic effects 



In the analysis, all secondary charged particles have been assumed to be pions, unless 
identified as muons, electrons or protons. However, as seen in section 3 for simulated 
events, there is a fraction of 39% and 23%, respectively of these positive and negative 
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particles, which are not pions. Here we want to study how the BEC could be changed by 
these misidentified tracks and by the use of the unlike-sign sample as a reference. These 
effects could manifest themselves in three distinct ways: 

• BEC for like-sign kaons and correlations for like-sign fermions (electrons, protons). 
The BEC for kaons were measured at LEP and showed characteristics very similar to 
those of pions. Fermion pairs, instead, could exhibit an anticorrelation effect. However, 
the number of like-sign kaon and fermion pairs is very small and their contribution is 
negligible. 

• Pairs of like-sign, but not identical particles in the numerator of eq. 5, for example 
K + 7i + pairs, have no BEC. This contribution could bias the value of A. 

• Unlike-sign pairs in the denominator of eq. 5 include contributions from K° and reso- 
nances such as p, u, as well as from electron-positron pairs from photon conversions. 
The latter could severely affect the distributions of R at small Q. 

Fig. 3 (a) shows the simulated Q distributions for reconstructed particles associated to 
generated primary pion pairs for like and unlike-charge distributions. One can notice that 
the unlike-sign pair distribution exhibits strong enhancements around Q ~ 0.35 GeV 
and Q ~ 0.7 GeV due to K° and p decays, respectively. For this reason the Q intervals 
0.3 ^ Q ^ 0.45 and 0.6 ^ Q ^ 0.825 GeV have been excluded from the analysis. The 
ratio like/unlike pairs, eq. 5, is also shown in fig. 3: as expected, no structure is observed 
at low Q. Adding kaons and protons we obtain the Q distribution shown in fig. 4 (a). 

W) 



0.2 0.4 0.6 0.8 1 1.2 1.4 " " .2 0.4 0.6 0.8 1 1.2 1.4 

Q (GeV) q (QeV) 

Fig. 3. (a), Q- distributions for Monte Carlo generated pion pairs; top to bottom: unlike, like, (+ 
+), ( — ) pairs, (b), the ratio like/unlike pairs, R(Q) of eq. 5, for a pure pion MC sample. 

The Q variable was calculated assigning the pion mass to all particles. No structure is 
visible, only a global shift towards higher values. Fig. 4 (b), shows the effect of adding 
all other particles including the "not recognized" ones in the sample: the first bin is now 
low, demonstrating that the denominator of eq. 5 contains a sizeable contribution at very 
low values of Q due to e + e~ pairs from photon conversions. For this reason the data at 
Q ^ 0.04 GeV have been excluded from the fit used to extract the BEC parameters. 

We observe that the unlike-sign pair distribution as reference sample has the essential 
property of reproducing faithfully the non-BEC distribution of like-sign pairs (the ratio 
is flat). However, it is dangerously affected by meson resonances and by electron-positron 
pairs from photon conversions. In particular the conversions give a major contribution to 
our systematic errors which will be estimated in section 6. 
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From this study we conclude that contaminations from particles other than electrons 
produce only a variation of the overall normalization of the distributions and no distortion 
of its shape. Therefore they do not affect the measurement of the radius of the emitting 
source while some effects could be induced on the chaoticity parameter. The observation, 

described in section 5, that the results obtained for the ( ) sample are very similar to 

those obtained for the (+ +) pairs, although the two samples are affected by different 
contaminations, demonstrates that these latter are not a critical issue. 



• 1 

0.8 
0.6 
0.4 
0.2 



0.2 0.4 0.6 0.8 



1.2 1.4 

Q (GeV) 



0.6 



0.2 





(b) 






_+ 








) 0.2 0.4 0.6 0.8 


1 1.2 1.4 



Q (GeV) 



Fig. 4. (a), R(Q) for a Monte Carlo sample of pions, kaons and protons, (b), R(Q) for a Monte 
Carlo sample of all particle types. 



5 Results 



This section presents the results on the chaoticity parameter A and the source radius R 
obtained following the Goldhaber, KP and (Q\\, Q±) parametrizations. 



5. 1 The Goldhaber parametrization 



The inclusive experimental correlation R(Q) as a function of Q is shown in fig. 5. The 
empty regions in the distribution correspond to the excluded intervals described above. 
Superimposed to the data is a fit of the form: 

R(Q) = N[l + Xexp(-R 2 G Q 2 )](l + aQ + bQ 2 ) (6) 

where N is a normalization constant and the second degree polynomial is a parametriza- 
tion of the shape of the long-range correlations outside the BEC region. The choice of the 
parametrization used to describe the long-range correlations inevitably affects the results 
of the BEC analysis and contributes to the systematic errors on A and Rq- The second 
degree polynomial gives the best x 2 /d.o.f. (compared to a quadratic or linear long-range 
form) of the fit and it has been often used in other experiments. Therefore, we shall use 
it in this paper and we shall discuss the use of other parametrizations in the section on 
systematic uncertainties. 
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Fig. 6 shows the BEC for (+ +) pairs (a) and ( — ) pairs (b) again with a fit of the form 
of eq. 6 superimposed. The long-range correlations are very different between positive 
and negative pairs: for positive pairs there is a steady increase of the correlation function 
for Q ^ 0.4 GeV, while for negative pairs the correlation function is almost flat at large 
Q. The long-range correlations for positive pairs also determine the behaviour of the 
like-sign pair sample. However, as opposed to the long-range correlations, the BEC are 

similar for the (++) and ( ) samples. Table 2 summarizes the results on A and R G for 

the like, (++) and ( ) samples. The values of A and R G of the three samples are in 

good agreement, in spite of the different shape of the long-range correlation region. This 
demonstrates that our parametrization (6) is robust and capable of describing correctly all 
three sets of data. The BEC parameter A is about 0.4 and R G is about 1 fm, independent 
of the particle charge. We notice the following: 

• The measured R distribution, in the region Q ^ 0.5 GeV (the region of the long-range 
correlation) differs from the Monte Carlo simulation, which produces a flat distribution 
(see fig. 4). 

• The best fit x 2 value is inconsistent with statistical errors alone. Most of the contribution 
to the x 2 comes from the region of long-range correlations (Q ^ 0.8 GeV) which is 
not fully accounted for by our empirical parametrization. Results from the NOMAD 
experiment on the production of the /o(980) and ^(1270) resonances in interactions 
have been published [19]. These resonances contaminate the region at Q > 0.8 GeV and 
contribute to the large x 2 of the fit. We verified that the exclusion of the region 0.9 < 
Q < 1.3 GeV improves the x 2 /d.o.f. from 1.75 to 1.5 and does not affect significantly 
the results. Moreover in the BEC region (Q < 0.2 GeV) the quality of the fit is always 
good. 



Pairs 


A 


Rg (fm) 


X 2 /d.o.f. 


like 


0.40 ± 0.03 


1.01 ±0.05 


90/52 


(+ +) 


0.38 ± 0.04 


1.03 ±0.07 


80/52 


(— ) 


0.43 ± 0.04 


0.96 ±0.06 


75/52 



Table 2 

Chaoticity A and Goldhaber radius R G . Errors are statistical only. 



5.2 The Kopylov-Podgoretskii parametrization 



The dependence of R on Q t and Qq is shown in fig. 7 (a) for ( ) pairs. A peak at Q t 

and Q ~ is visible with a width of 0.2 GeV in both variables. 

To fit the two dimensional structure of BEC it is convenient to use a one dimensional 
Qt parametrization derived from eq. 3 (see [6], [17]) by restricting the allowed energy 
difference to Q ^ Q ma x- Then, under the hypothesis that (Q m axT~) 2 <C 1: 



C{Q t ) = N KP 1 + \ K p[2Ji(R K pQt)/(RKpQt)} 2 (1 + aQ t + bQ 2 t ) 



(7) 
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Q_4 I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 

0.2 0.4 0.6 0.8 1 1.2 1.4 

Q (GeV) 



Fig. 5. R as a function of Q in the like-sign pair sample. Superimposed is a fit obtained using 
the Goldhaber parametrization (eq. 6). 




Fig. 6. R as a function of Q in the (++) (a) and ( ) (b) pair sample. Superimposed is a fit 

following the Goldhaber parametrization (eq. 6). 

Here too a polynomial form is used to parametrize the long-range correlation. The Q t 
distribution for Qq ^ 0.2 GeV is shown in fig. 7 (b). Notice that also in this case we 
remove from the fit the two regions where the presence of K° and resonance contributions 
affects the Q t variable. 

The result of a fit using the parametrization given in eq. 7 yields Xkp = 0.29 ± 0.06 and 
Rkp = 2.07 ± 0.04 fm, in agreement with the expected relation Rrp ~ 2Rg- Again the 
results obtained using the (++) and ( ) samples are consistent with each other. 
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Pairs 




Rkp (fm) 


X 2 /d.o.f. 


like 


0.29 ±0.06 


2.07 ±0.04 


51/52 


(+ +) 


0.28 ± 0.04 


2.13 ±0.04 


56/52 


( — ) 


0.32 ±0.06 


2.01 ±0.04 


38/52 



Table 3 

Chaoticity parameter Xkp and radius Rkp- Errors are statistical only. 




Fig. 7. (a), R as a function of the KP variables Qt and Qq for ( ) pairs, (b), R as a function 

of Qt for Qo ^ 0.2 GeV. Superimposed is a fit using eq. 7. 



5.3 Longitudinal and transverse shapes 



A possible deviation of the pion emitting source from a spherical shape in its rest frame can 
be investigated in the LCMS system using the variables Q± and Qn . The BEC behaviour is 
studied separately for the two variables by requiring Q» ^ 0.2 GeV for the Q± distribution 
and, conversely, Q± ^ 0.2 GeV for the Q\\ distribution. These distributions are shown in 
Fig. 8. Again the regions where the presence of K 0, s and resonances affects Q» and Q± 
have been removed from the fit. The fit is performed with a parametrization as in eq. 
4, multiplied by a second degree polynomial to reproduce the long-range correlations. 
The fitted values for the BEC parameters are shown in tables 4 and 5 together with the 

results obtained from (++) and ( ) pairs separately. Our measurements confirm the 

LEP results [22] that in the LCMS reference frame the longitudinal size of the pion source 
is 30-40% larger than the transverse one. 
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Pairs 


R± (fm) 


X 2 /d.o.f. 


like 


0.98 ±0.10 


71/52 


(+ +) 


1.04 ±0.12 


63/52 


(— ) 


0.81 ±0.15 


50/52 



Table 4 

The LCMS variable R± for Qy ^ 0.2 GeV. Errors are statistical only. 



Pairs 


R\\ (M 


X 2 /d.o.f. 


like 


1.32 ±0.14 


54/52 


(+ ±) 


1.39 ±0.24 


50/52 


(— ) 


1.15 ±0.12 


64/52 



Table 5 

The LCMS variable Ru for Q± ^ 0.2 GeV. Errors are statistical only. 
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Fig. 8. (a), R vs Q± for Q\\ ^ 0.2 GeV. (b), R vs Q\\ for Q± ^ 0.2 GeV . Superimposed is a fit 
using eq. 4 multiplied by a second degree polynomial. 



5.4 BEC dependence on the rapidity of the pair 



Deep inelastic CC neutrino interactions involve a d quark in the target nucleon leaving 
as spectators the remaining quarks. We naively expect, therefore, two distinct sources 
of secondary hadrons: the single struck d quark and the spectators. At high energy the 
two contributions should be fairly well separated in the cm. frame of the hadronic jet 
as the particles coming from the fragmentation of the struck quark should have positive 
rapidities while those produced by the spectator quarks should have negative rapidities. 
To investigate possible differences between the two pion sources we studied the BEC 
distributions, using the Goldhaber parameter Q, for pairs of particles of equal rapidity 
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sign and also for pairs of particles of opposite rapidity sign in the rest frame of the hadronic 
jet. The data are shown in fig. 9 and 10. Table 6 summarizes the results obtained on A 



Rapidity- 


A 


R G (fm) 


X 2 /d.o.f. 


Inclusive 


0.40 ± 0.03 


1.01 ±0.05 


90/52 


Positive 


0.47 ± 0.04 


0.98 ±0.07 


84/52 


Negative 


0.42 ± 0.09 


1.03 ±0.17 


53/52 


Opposite 


0.37 ±0.06 


0.98 ±0.08 


52/52 



Table 6 

Rapidity dependence of the chaoticity A and the Goldhaber radius Rq- Errors are statistical only. 

and Rq for the various rapidity configurations. The source radius Rq shows no differences 
for particles emitted at different rapidities, demonstrating that the typical hadronization 
scale is much longer than the interaction radius, resulting in a unique hadron source, 
independent of the detail of the quark interactions. 



R 

1.2 
1.1 



0.9 
0.8 
0.7 
0.6 
0.5 
0.4 

0.2 0.4 0.6 0.8 1 1.2 1.4 

Q (GeV) 

Fig. 9. R vs Q for positive rapidity sign pairs. Superimposed is a fit using eq. 6. 



5.5 BEC dependence on event charged multiplicity and hadronic energy W 

The large number of events collected by NOMAD allows the study of BEC effects in dif- 
ferent final state configurations. In particular it is interesting to verify the observations in 
hadronic [26] and e + e~ interactions [20] that the Goldhaber radius increases with the event 
charged multiplicity N c h- Fig. 11 (a) shows the Goldhaber radius Rq and the chaoticity 
parameter A for seven different N c h values. We see that there is here an indication for a 
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Fig. 10. R vs Q for negative (a) and opposite (b) rapidity sign pairs. Superimposed are fits using 
eq. 6. 

decrease of Rc with N c h- One should notice that the rise of the emission radius Rq with 
N c h at LEP is only visible at very high multiplicities (N c h > 10) which are not accessible 
to this experiment. The chaoticity parameter A appears to increase with N c h- 




Fig. 11. Dependence of the chaoticity parameter A and Goldhaber radius Rg on the event charged 
multiplicity N c h (a) and on the hadronic energy W (b). 



A similar conclusion can be drawn when studying the BEC effects as a function of the 
hadronic energy W: fig. 11 (b) presents Rq and A for six W intervals. Rq decreases with 
W whereas A increases. 
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6 Systematic errors 

We focus the discussion of systematic errors to the inclusive BEC study using the Gold- 
haber parametrization. 

Coulomb interactions between particles which affect like-sign and unlike-sign pairs in 
opposite ways, can alter the correlations. This effect changes the two pion cross section 
by the Gamow factor [16], a significant correction only at very small values of Q. We 
checked that the effect enhances both A and Rg by only a few percent therefore we 
decided not to apply it. 

We identify three sources of systematic errors in our results: 

• The uncertainty in the background contribution underneath the BEC peak at Q < 
0.2 GeV. In the NOMAD experiment BEC effects could be affected by an insufficient 
rejection of e + e~ pairs from photon conversions contaminating the reference sample. 

• The cuts applied. 

• The track reconstruction efficiency We verified that this effect does not produce sizeable 
effects on the Goldhaber parameters, decreasing the track reconstruction efficiency by 
10%. This was done by removing tracks at random from the sample before calculating 
the BEC. 

6. 1 Systematic errors from the e + e~ background 

As already noted in section 4.1, BEC could be altered by the presence of background of 
e + e~ pairs from photon conversions in the unlike-sign sample used as a reference. Most 
of this background is at low Q and could seriously affect the results. To reduce this effect 
the data at Q ^ 0.04 GeV in all the previous correlation distributions have been excluded 
from the fit used to extract the BEC parameters. The systematic uncertainty from this 
cut is then estimated by enlarging the data excluded from the fit to Q ^ 0.06 GeV 
(i.e. the first bin of fig. 5). The results obtained for A and Rg are shown in table 7. By 
comparing table 7 to table 2 we see that our results are insensitive to a variation of the 
lowest accepted Q bin. A larger sensitivity of A and Rq is found when varying the track 
quality cuts (see next paragraph). 



Pairs 


A 


Rg (fin) 


X 2 /d.o.f. 


like 


0.40 ± 0.03 


1.02 ±0.05 


89/51 


(+ +) 


0.39 ± 0.03 


1.04 ±0.06 


80/51 


( — ) 


0.44 ±0.05 


0.98 ± 0.08 


84/51 



Table 7 

Chaoticity parameter A and Goldhaber radius Rq obtained from a fit to the R distribution where 
the data at Q ^ 0.06 GeV have been removed. 
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6.2 Systematic errors due to the selection cuts 

We checked the stability of our results by varying the following track selection parameters: 

• The difference Az between the first point of the reconstructed track and the primary 
vertex position along z: Az $C 15 cm; 

• the minimum track momentum: \p\ ^ 100 MeV/c; 

• the maximum acceptable momentum uncertainty ^ 6%. 

We note that these selection parameters affect differently pions and electrons and, there- 
fore, they could change the fraction of conversions included in the data. 

To estimate the effect of varying these cuts on the fitted parameters A and Rc each cut 
was modified and the relevant BEC distribution was again fitted. The results are shown 
in table 8. The interval chosen for Az corresponds to the thickness of one DC chamber. 
The variations of the momentum and cut positions in table 8 reflect the uncertainty 
in their choice for an optimum separation of the electron and the pion populations. A 
similar procedure was adopted also for ( ) and (+ +) pairs separately. The largest 



cut 


A 


R G (fm) 


all cuts 


0.40 ±0.03 


1.01 ±0.05 


Az 10 cm 


0.40 ±0.04 


1.10 ±0.09 


Az ^ 20 cm 


0.37 ±0.03 


0.95 ±0.04 


\p\ ^ 0/0 


0.35 ±0.03 


0.99 ±0.04 


|p| ^ 150 MeV/c 


0.40 ±0.03 


1.01 ±0.05 



Table 8 

Chaoticity parameter A and Goldhaber radius Rq obtained for like pairs and for different cut 
configurations. 

effect is the one induced on Rq by changes in the Az cut and on A by changes in the ^0- 
cut. The systematic uncertainty due to the Az cut can be also estimated by extrapolating 
the parameters to Az = 0. The variations for Rc and A amount to 15% and to 12% 
respectively. 

6.3 Effect of the long-range correlation parametrization 

The numerical values of Rq and A depend on the parametrization used to describe the 
long-range correlations. In the literature linear, quadratic and polynomial forms have 
been used. This ambiguity must be taken into account when comparing results from 
different experiments. Throughout this paper the long range correlations, necessary to 
describe effects other than BEC, have been described by a second degree polynomial 
form (1 + aQ + bQ 2 ). Different parametrizations are possible and have been used in other 
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experiments: i.e. a linear form: (1 + aQ) or a quadratic form: (1 + bQ 2 ). Table 9 shows 
the results of the three different choices for the long range parametrizations when the like 
distribution is analyzed. We see that the fit worsens when using the linear parametrization 
and there are also significant differences in the fit values for A and Rq, while the quadratic 
parametrization reproduces almost exactly the results of the polynomial. We conclude 
that the linear approximation is inadequate to be used in the analysis of these data. 
For completeness in section 8 our results obtained with different parametrizations of the 
long range effects will be compared with the data of other experiments using similar 
parametrizations. 



Long range parametrization 


A 


Rg (fin) 


X 2 /d.o.f. 


Polynomial 


0.40 ± 0.03 


1.01 ±0.05 


90/52 


Linear 


0.54 ± 0.02 


0.86 ± 0.03 


123/52 


Quadratic 


0.43 ± 0.02 


0.95 ± 0.03 


93/52 



Table 9 

Chaoticity parameter A and Goldhaber radius Rq for three choices of the long range correlation 
parametrization in the NOMAD analysis. Errors are only statistical. 



7 Final results 

Table 10 summarizes our final results on A and Rq including also the systematic errors 
from variations of the cuts discussed in the previous section (added in quadrature). 



Pairs 


A 


Rg (fm) 


like 

(+ +) 
(— ) 


0.40 ± 0.03±g;g| 
0.38 ± 0.04±g;{£ 
0.43 ± 0.04lg$ 


i.oi ± o.o5±8:gg 

1.03 ± 0.07lg;g? 
0.96 ± 0.06±g;^ 



Table 10 

Chaoticity A and Goldhaber radius Rq. The first error is statistical, the second one is systematic. 



8 Comparison with the results of other experiments 

Fig. 12 and table 11 display a compilation of some measurements of A and Rq in the 7T7T 
channel in high statistics lepton-induced reactions: neutrino interactions [6], muon DIS 
[18], electron-proton DIS [21], [27], e+e~ collisions [23], [24], [20], [25]. These experiments 
were performed at different energies; they have different selection criteria and biases and 
also different parametrizations for the long-range correlation (see table 11). 
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Our results agree within errors with those of the combined analysis the BEBC and of the 
Fermilab neutrino data [6]. 

The results on A (fig. 12 (a)) shows that there are two groups of experiments which are 
consistent within each group, but not between them. The first group clusters around 
A fa 0.5 and the other around A ps 1. The parametrization of long-range correlations does 
not seem to be the origin of this discrepancy. Since A is sensitive to the purity of the pion 
sample the origin of this difference could be the different pion identification criteria of the 
experiments. 

The results on R G are shown in fig. 12 (b): it appears that the value of R G computed 
with a linear model is systematically lower than the one computed with a quadratic or 
polynomial form. The two groups of data are rather well consistent within each other: the 
"linear" group cluster at R G ~ 0.6 fm and the "quadratic-polynomial" group at R G ~ 0.9 
fm. 



Experiment 


Fit type 


(GeV) 


(Q 2 ) 

(GeV 2 ) 


A 


Rg 
(fm) 


BEBC-Fermilab 


quadratic 


10 


10 


0.61 ±0.04 ±0.15 


0.80 ±0.04 ±0.16 


EMC 


quadratic 


23 


50 


1.08 ±0.1 


0.84 ±0.03 


DELPHI 


linear 


91 




1.06 ±0.05 ±0.16 


0.49 ±0.01 ±0.05 


ALEPH 


linear 


91 




0.51 ±0.04 ±0.11 


0.65 ±0.04 ±0.16 


ZEUS (DIS) 


linear 


300 


400 


0.431 ± 0.012l°;^ 2 


0.671 ± 0.016lao32 


HI (DIS) 


linear 


300 


40 


0.52 ± 0.031&21 


0.68 ± 0.04lg;^ 


L3 


linear 


189 




0.48 ± 0.05 ± 0.07 


0.71 ± 0.04 ± 0.05 


OPAL 


polynomial 


91 




0.672 ±0.013 ±0.024 


0.955 ±0.012 ±0.015 


NOMAD 


polynomial 


8 


10 


0.40 ± 0.03±g;ge 


1.01 ± 0.0518;" 


NOMAD 


quadratic 


8 


10 


0.43 ± 0.02 


0.95 ±0.03 


NOMAD 


linear 


8 


10 


0.54 ±0.02 


0.86 ±0.03 



Table 11 



Summary of results published in previous experiments. 



9 Conclusions 

The NOMAD experiment has measured BEC in charged- current neutrino interactions 
using different parametrizations for this effect. The general picture emerging from the 
data is that the size and the chaoticity of the pion source are about 1 fm and about 
0.4 respectively, quite independent of the final state rapidity sign of the emitted pions. A 
difference of about 40% is found between the longitudinal and transverse size of the source. 
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Fig. 12. Compilation of results obtained by various experiments for the chaoticity parameter A 
(a) and the Goldhaber radius Rq (b). 

We observe a decrease of the Goldhaber radius as a function on the charged multiplicity 
and of the hadronic energy of the event. A comparison of our results with those of other 
experiments studying other processes than neutrino interactions shows a fair agreement, 
demonstrating that the final state hadronization processes have universal features with 
little dependence on the type or energy of the interacting particles. 
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